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Abstract 

Nonlinear evolution equations of the fourth order and its partial 
cases are derived for describing nonlinear pressure waves in a mix- 
ture liquid and gas bubbles. Influence of viscosity and heat transfer 
is taken into account. Exact solutions of nonlinear evolution equa- 
tion of the fourth order are found by means of the simplest equation 
method. Properties of nonlinear waves in a liquid with gas bubbles 
are discussed. 

1 Introduction 

It is well known that a mixture of liquid and gas bubbles can be considered as 
an example of a nonlinear medium. Van Wijngarden [1] was first who derived 
that the famous Korteweg - de Vries equation [2H1] can be used for describing 
nonlinear waves in liquid with gas bubbles. Taking into account the viscosity 
of liquid, Nakoryakov and et al. [S] obtained the Burgers equation [HIE] and 
the Korteweg - de Vries - Burgers equation to describe the pressure waves in 
bubbly liquids. The dynamic propagation of acoustic waves in a half - space 
filled with a viscous, bubbly liquid under van Wijngaarden linear theory was 
considered in the recent work [8J. Oganyan in [H] investigated a quasi - 
isothermal wave propagation in a gas - liquid mixture, in which he take into 
consideration the thermodynamic behavior of the near - isothermal gas in 
the bubbles. 
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Many authors studied nonlinear processes in a liquid with gas bubbles 
using the numerical methods. For an example Nigmatullin and Khabeev 
[TQl |TT] considered the heat transfer between an gas bubble and a liquid 
and studied the structure of shock waves in a liquid with gas bubbles with 
consideration for the heat transfer. 

The purpose of this work is to obtain nonlinear evolution equations for 
describing the pressure waves of in a liquid with gas bubbles taking into 
account the viscosity of liquid and the heat transfer on boundary a bubble 
and a liquid. 

2 System of equations for describing waves 
in a liquid with gas bubbles with consider- 
ation for heat transfer and viscosity 

Suppose that a mixture of a liquid and gas bubbles is homogeneous medium 
|T5] . In this case for description of this mixture we use the averaged 
temperature, velocity, density and pressure. We also assume that the gas 
bubbles has the same size and the amount of bubbles in the mass unit is 
constant. We take the processes of the heat transfer and viscosity on the 
boundary of bubble and liquid into account. We do not consider the processes 
of formation, destruction, and conglutination for the gas bubbles. 

We have that the volume and the mass of gas in the unit of the mass 
mixture can be written as 

V = ^nR 3 N, X = V Pg , 

where R = R(x,t) is bubble radius, N is number of bubbles in mass unit, 
p g = p g (x, t) is the gas density. Here and later we believe that the subscript g 
corresponds to the gas phase and subscript / corresponds to the liquid phase. 

We consider the long wavelength perturbations in a mixture of the liquid 
and the gas bubbles assuming that characteristic length of waves of pertur- 
bation more than distance between bubbles. We also assume, that distance 
between bubbles much more than the averaged radius of a bubble. 

We describe dynamics of a bubble using the Rayleigh - Lamb equation. 
We also take the equation of energy for a bubble and the state equation for 
the gas in a bubble into account. The system of equation for the description 
of the gas bubble takes the form [TJHIE] 

Pl (RR tt + -it!, 2 + ^Lll\ =Pg-P, (1) 
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P a ,» + ^, + 3 ^rV g -r,) = o, (2) 



Ts= np 1 fR\\ (3) 

where P(x, t) is a pressure of a gas-liquid mixture, P g is a gas pressure in 
a bubble, T 5 and 7] are temperatures of liquid and gas accordingly, Xg is a 
coefficient of the gas thermal conduction, Nu is the Nusselt number, n is a 
polytropic exponent, z/ is the viscosity of a liquid. 

The expression for the density of a mixture can be presented in the form 



Considering the small deviation of the bubble radius in comparison with 
the averaged radius of bubble, we have 

R(x, t) — Rq + T](x, t), Rq = const, ||r7||<<i?0) R(x,Q) = Rq. (5) 

Assume that the liquid temperature is constant and equal to the initial 
value 

Ti = T \ t=0 = T , T = const. 
At the initial moment, we also have 

t = 0: P = P g = To, = const, V = V = -nR 3 N. 

Substituting P g and T g from Eqs. ([1]) and (jHJ) into Eq. (j2J) and taking 
relation (J5]) into account we have the pressure dependence of a mixture on 
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the radius perturbation in the form 
P-P + £-P + s-Pvt + xP t + xpi(i? + v)Vttt + (3n + A)xpir] t r] tt + 



Pi(3R% + 4z/x) Pl (6R 2 Q - 4z/x) p ; (8z/x(3n - 1) + 9R 2 ) 

'It - 

(6) 



3R ^ + SR* mtt + 6^ % + 



JLvm 2P 3P 2 n 



x — 



3 X Nu(n - 1)T„ ' 

From Eq. (J3J) we also have the dependence p on 77 by means of formula flSJ) 
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p = po - fll] + ^r] , p 



l-X + ^p;' 

(7) 

3pjVo _ 6p^o(2p^ -l + X) 

^ i?o(i-^ + w ^ ^(l-x+p^o) 3 ' 

We use the system of equations for description of the motion of a gas- 
liquid mixture flow in the form 

dp d(pu) (du du\ dP 

dt dx ^ \ dt U dx J dx 

where u = u(x, t) is a velocity of a flow of a gas-liquid mixture. 

Eq.(jSJ) together with Eqs.flU]) and (j7j) can be applied for describing nonlin- 
ear waves in a gas-liquid medium. This system of equations can be written 



4 



as 



, Po 2/ii 
Vt + uVx + r\u x u x r]7]t = 0, 



-— (u t + UU X ) + 7]U t - -P x = 0, 



p + ~s~P + s-Pm + xPt = -m(Ro + v)vm- 

Pl (?>Rl + AvH) Pl (GRl-4vx) 
-(3n + ±)Kpir} t r)tt % ^2 Wht~ 



(9) 



p z (8^x(3n - 1) + 9-Rg) 2 4z/p, 2P 3P ( 



Consider the linear case of the system of equations Assuming, that 
pressure in a mixture is proportional to perturbation radius, we obtain the 
linear wave equation for the radius perturbations from Eqs.Q 



Vtt = eg Vxx, c = J (10) 
Let us introduce the following dimensionless variables 

t = — t', X = lx', U = C u' ', 7] = R 7] , P = P P' + P , 

Co 

where I is the characteristic length of wave. 

Using the dimensionless variables the system of equations can be 
reduced to the following (the primes are omitted) 

Po . . 2/iiPo n 
Vt ^~Ux + ur] x + r]u x r}T] t = 0, 

(u t + uu x ) + r]u t - -P x = 0, 



uRn 3 

(11) 



P + >c x P t + rjP + 3nxi rj t P = -jrjttt - JVVttt - (3n + 4)777^- 

i-L 3 
2 + 2' 



-(A + 02) Vtt - (2/3 2 - A) 7777^ - I — — + -fa ) vl 



-\r)t - 3r] + 3r] 2 , 
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where the parameters are determined by formulae 

4^co _ Aux Pl cl _ Pl c\ R\ 



( 12 ) 

7= -f^' X1= — 

Numerical and analytical task by means of the system of equations ( 11 ip 
is a difficult problem. 

3 Basic nonlinear evolution equation for de- 
scribing nonlinear waves in liquid with bub- 
bles 

We can not find the exact solutions of the system of nonlinear differential 
equations (ITT]) . To account for the slow variation of the waveform, we intro- 
duce a scale transformation of independent variables [17] 

£ = e m (x-t), T = e m+1 t, m>0, £<1, (13) 



dx dt ' <): d£ 



Substituting (TTBj) into (ITT]) and dividing on e m in first two equations we 
have the following system of equations 



er] T —r)£ — w § + i]U£ + ur]£ - e i]r] T H r]T]£ = (14) 



£ J ^ + ^ + ^ - ^ - - 3 Pf = ° (15) 
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P + £ m+1 xiP T - e m K l P i + r/P + £ m+1 3nxi?7 r P - e m 3nxi^P = 

-e^T^W + £ 3m+2 3 7 77^ - e 3m+1 3 7 r7^ + £ 3m 7 W^- 

-e 3m+3 (3n + 4) 7 r/ r r/ TT + e 3m+2 (3n + 4) 7 r/^ TT + 

+e 3m+2 2(3n + 4)7^77^ - e 3m+1 2(3n + 4)7^77^- 

-£ 3m+1 (3n + 4) 7 ^ a + £ 3m (3n + 4)777^ - e 2m+2 (ft + ft) Vrr + 

+ e 2m+1 2 (ft + ft) ^ - e 2m (ft + ft) m - £ 2m+2 (2/3 2 - ft) r/r/ rr - 

+ e 2m+1 2 (2/3 2 - ft) 7777^ - s 2m (2fa - fa) wit ~ 

~e 2m+2 + f&) r£ + e 2 ^2 + ffc) ^- 

/ Sr? — 1 3 \ 

~ £2m ( ~^2~ 01 + 2~N ^ ~ £m+1Xr lr + e m X^-3 V + 3 V 2 . 

Now we assume that the variables u, rj and P can be represented asymp- 
totically as series in powers of e about an equilibrium state 

u = eu\ + e 2 u 2 + • • • , r] = er]i+ e 2 T] 2 + P — eP\ + £ 2 P 2 + . . . (17) 

Substituting ffTTl) in fTH|) - f[TB]) and equating expressions at e, we have the 
equations 

Solving this system of equation, we have 

Vx(Z,t) = — ^-Ui(£,t) + ^(t), P 1 (e,r) = -377 1 (e,r), (18) 

where ?/>(t) is an arbitrary function. We suppose below that ip(r) = 0. 

Substituting fJTTji in (fl4|) - (fr6|) and equating factors at e 2 to zero, we have 
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the system of equations 

Po %UiRo 

Vir - ^~U2( + Vi u U + UiTjit H r]ii] U = 0, 

fiKo fj, 



Po Po Po 1 P n 



pRo pRo pRo 3 

(19) 

P 2 - e™- 1 ^ + TfcPi - e m 3nx 1 »^P 1 = e^-^K + e 3m 7»7i»te«+ 
+e 3m (3n + 4)^^ - e 2 ™" 1 (A + &) Vm - e 2m {2fc - (3,) VlVm - 
3n- 1 



—e 2m 



+ -ft) ^ + e^A ^ - 3^ + 377?. 



Taking Eqs. flTB"]) and ( TT9|) into account we have the equation for pressure 
Pi in the form 



p 4- n P P _l_ P 2m-1 ^1 p 2m ^2 ~ Pi p p 

f Xr +«A^+£ — - — f im - e — — — n% - 



i_m ^"^1 ( td jd \ i _3m— 1 7 p -3m 7 p p 

+ £ "Tr - lAPcJf + £ 7T Pl£5« - £ 777 A -^i^ - 



'^ T ' ii 18 '"^ (20) 



^3m (gw + 5)7 3m (3n + 4) 7 



3fiR() 3/iii?o 2 

a = + -. 

Po P 3 

Assuming m = 1 we obtain the Burgers equation from Eq. f[2"Uj) when 
e -»■ 

P lT Q - y) (21) 

Many solutions of the Burgers equation are well known [3J[T5J[ni]. By 
means of the Cole - Hopf transformation [201ET] the Burgers equation can 
be transformed to the linear heat equation. The Cauchy problem for this 
equation can be solved [3]. 
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In the case m — |, (| — ^) ~ e 2 we have the famous Korteweg - de 
Vries equation [2J at e — )■ 

P lT + a P l P li + P im = 0. (22) 

This equation has the soliton solutions [22] . The Cauchy problem for the 
Korteweg - de Vries equation can be solved by the inverse scattering trans- 
form [23] • There are many methods to find exact solutions of the Korteweg 
- de Vries equation [2H[25]. Unfortunately some of these approaches lead to 
the erroneous solutions [26H28] . 

Assuming m = |, (| — ^) ^ e we have the Korteweg - de Vries - 
Burgers equation at e — » 

P lr + a P lPl( + P lm = g - f ) fl« (23) 

In the case m = |, /3i + /3 2 — 6 e^, (~ — ^) — e ^ we have the famous 
Kuramoto - Sivashinsky equation [29| l30] when e — > 

P lr + a PxP lf + P 15C€ = P 1K + 1 P aKK . (24) 

The cauchy problem for equation f[2~4"|) can not be solved by the inverse 
scattering transform but this equation has some exact solutions. The Ku- 
ramoto - Sivashinsky equation was studied in many papers (see, for the ex- 
ample [aMH]). 

From Eq. fl20p one can find some other nonlinear evolution equations for 
describing waves in a mixture of liquid and gas bubbles. 



4 Fourth order equation for describing wave 
processes in liquid with gas bubbles 

Assuming m = ~, A = 6e 2 ^ 3 \' , A' ~ 1, K\ = 2e 2 ^x 1 , x x ~ 1, f3 1 = 6s 1 / 3 (3^, 
P[ ~ 1, 2 = 6e^ 3 (3' 2 , /?2 ~ 1, S > in Eq.fl20]) and usine; transformations 

Pi = - P[,a = e a we have nonlinear evolution equation of the fourth order 
in the form 



pi + « pi?* + (/?; + a) ^ - (^^) ^« - 

0[ + = (a' - xi) p; K + s (p; f p; 



3n - 2 



(25) 



rip ^ p' P > (3n + 5) 7 , , (3n + 4) 7 , 2 
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Coefficient 7/6 at the fourth derivative in (|25p can be presented in the 
form 

r = 1 = 4 ^ C o^ (26) 

6 9 Xg T Nu{n-l)l 3 ' 1 ; 



We can see that the value T in (126]) depends on the initial radius of 
a bubble, sound velocity, coefficient of the heat transfer and the Nusselt 
mumber. The heat capacity of liquid is much more than the the heat capacity 
of the gas in bubble. As consequence the temperature of liquid is not changed 
but the gas in bubble can be heated and cooled. In the case of the isothermal 
processes we have Nu — > 00 and Eq. (l25|) becomes the third order equation. 

Let us look for the exact solutions of equation (1251) . Using the variables 

c - 1 C ' r= t r ' P > = ZM±M V 

? Q((3'i + fc)^ 216(# + #)< ' 1 2#-/9i ■ 

we obtain the equation (the primes are omitted) 



V T + Oil V V £ + Vtf£ - (yVtf)^ - /3 V£V££ = 
= <TV I £ + Oi {VV^^ + V aa - 7 2 (UU«{){ - 73 {V(V^) ( 

where parameters of Eq. (j2"7|) can be written in the form 



(27) 
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a'l 2 (3n - l)(3'i + 3/3 2 7 (V - xQ 

ai 12(/3; + ^) 2 (2^-/3i)' P 2/3;-/3; ' ° 6(&+&) 

n^il (3'i + (3' 2 (3n + 4)(f3' 1 + f3' 2 ) 

°"l = TTTTTTv o'wn' , ^7T ; 72 = TTTy TT-, 73 



2(2& + " 2/?2 - ' J 2^-/3; 

Using traveling wave ansatz t>(£, r) = y(z), z = £ — C r and integrating 
the equation on z, we have 

Ci - C y + y y + y« - yy« - ^Vz ~ a y* ~ a i yy z - ^ 

-Vzzz + 72 + 73 UzVzz = 0. 

To find the solution of (|2"H|) we use the simplest equation method [IHH5I] 

y(z) = do + w, w = iw(z), (29) 
where u>(z) is the solution of the linear equation 

w zz = Aw z + Bw. (30) 
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Substituting the transformation and taking into account Eq. ([3D]) we 
obtain the following solution of the equation 



y(z) = a + C 3 e kz + C 2 e~ kz , (31) 

(/%3 + 72) + a/ (3 2 (73 + 72) (73 - 3 72) + 8 73(73 + 72) (/3 + 2 at 73) 
£ — _v 

473 (73 + 72) 

where C 2 and C3 are arbitrary constants, parameters ao, A, £?, «i Co and C\ 
are determined by formulae 

_ 2 g (72 + 73 - 1) - 4 a 73(73 + 72) - 4 q-j 73 - (3 2 
a °~ 273/3 + 4a l7 3 2 -72/3 2 

A= JL B _ -72/3 2 + 27 3 /3 + 4a 1 7 3 2 
273 ' 4 (73 + 72) 7! 

_ 4 7 2 /3 + 7 2 /3 3 - 4 72 /3 2 73 + 8 a x 7 3 3 - 4 ax 7 2 /3 72 
1 4 7 f ( 73 + 72) 

C = - 37 1 ■ r (4ai7 3 3 - 4(x 1 73 2 - 2(Ti73 2 /3 + 2(T 1 /372 73 + 72/3 2 + 

471(73 + 72) 

+2 a 7 3 2 /3 72 - -4 7 |a - 4 7 3 2 72 a + 4 7 2 /3 + 2 /3 72 73 - 2 /3 2 7 3 - 2 72 /3 2 7s+ 
+2 73 a/3 72 2 - 273/3 + /3 3 72 - 72 2 /3 2 ) , 

= Qfl- 2 73) (/3 2 -2/37 2 + 2/3- 273/3 + 472a73 + 4ai 73 + 4a7 3 2 ) 
1 87I (73 + 72) 

For finding the periodic solutions of Eq. (|28|) we use the expression in the 
form 

y(z) = ao + cti w + a 2 w 2 + a 3 w 3 , w = w(z) (32) 

where ^(2;) also satisfies linear ordinary differential equation (130]) . 

Substituting Eq. (l32|) into Eq. (f28|) and equating the different expressions 
at w (z) to zero, we have the system of the algebraic equations/, Solving this 
system of equations we obtain the following values of parameters 



1 Ot\ 

a 2 = a 1 = 0, a = 1, 72 = 1, 73 = -r, B = 



3' "~ 6(A-1) 
-2A 2 + A- ai + A 3 -2 A 2 + A- a x + A 3 

a = A=l ' ai = A=l 

10A-A ai 
P = ^ , 6=0, C = ai, 61 = — 



(33) 
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and solution of equation ( 12 8p in the form 



y(z) = 1 + a 3 



C 3 exp [ z — + . lN + C 2 exp z — 



2 6(A- 1) J z * \ 2 6(A-1 
p = 3(A-l)(3A 3 -3A 2 -2ai), 0, A ^ 1, a 3 ^ 

(34) 



In the case of 



(-3A + 3A 2 + ^) (-3A + 3A 2 - yp) 

6(A-1) ' 6(A-1) < 



solution (1341 is finite function. 

Consider case A = 0. In this case the parameters of equation and coeffi- 
cients of expansion take the form 

1 «i 

a 1= a 2 = 0, a = 1, 72 = 1, 73 = B = — 

3 o 

4 ax 
a = -ai, <Ti = «i, /? = --, C = 0, C = ai, Ci = y. 



At ot\ > we have the following solution of equation (|28|) 

y(z) = 1 + a 3 (c 3 + C 2 e~^^j ,a 3 ^0 (35) 

In the case of ol\ < the solution of equation (|28p can be written as 

V(z) = l + a 3 (a Bin {^a} + C 2 cos {^S}) 3 , «3 * 0. 

(36) 

Solutions (135|) . ( 136|) of Eq. (l28|) have two arbitrary constants C 2 and C3. 

Solutions (13"4"j) . ( 155|) have physical sense. The first of them is limited at 
the values of parameters at conditions of the above mentioned. The solution 
( 136]) is the periodical solution. The velocity of the wave described ( IMj) . (I35|) . 
(I36|) . is equal to parameter characterizing nonlinear transfer. 

Dependence of a solution (l36|) of Eq. ([28]) with respect to z at values of 
parameters a 1 = —1, C 2 = 0.04, C 3 = 0.05 is given on Fig. 1. 
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Figure 1: Periodic solution f[3"6"j) of Eq. ( 12"5|) 

5 Conclusion 

We have studied the propagation of nonlinear waves in a mixture of liq- 
uid and gas bubbles taking into consideration the influence of heat transfer 
and viscosity. We have applied different scales of time and coordinate and 
the method of perturbation to obtain nonlinear evolution equations for de- 
scribing pressure waves. As a result we have found the nonlinear evolution 
equation f l2"Uj) of the fourth order for describing the pressure waves in a liq- 
uid with bubbles. This evolution equation allow us to take into account the 
influence of the viscosity and the heat transfer on the pressure waves. It is 
likely this nonlinear differential equation is new and generalizes many well - 
known nonlinear equations as the Burgers equation, the Korteweg - de Vries 
equation, the Korteweg - de Vries - Burgers equation and the Kuramoto 
- Sivashinsky equation. Generally this equation is not integrable but this 
equation has some exact solutions that have been found in this paper. 
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